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1. Introduction 

Ternary fission involving the emission of a-particle was first observed |1[] more than fifty 
years ago. Emission of a-particles in the spontaneous fission of 252Qf Yia,s also a long 
history of investigation experimentally and theoretically p, ^ . 

A renewed interest in the spontaneous fission of ^^^Cf arose in connection with 
modern experimental techniques, ^ — ^ — ^y and x — 7 — 7 triple coincidence, 
(Gammasphere with 110 Compton suppressed Ge detectors), which allow the fine 
resolution of the mass, charge and angular momentum content of the fragments [§]. Very 
recently direct experimental evidence was presented for the cold (neutronless) ternary 
spontaneous fission of 252Qf j^ -^vhjch the third particle is an a-particle 0, or ^°Be 0. 
This confirms that a large variety of nuclear large-amplitude collective motions such 
as bimodal fission [§, cold binary fission 0, ||, |10], |Tl], 0, heavy cluster radioactivity 
||1^ , p!^ , and inverse processes, such as subbarier fusion ||T^, could belong to the general 
phenomenon of cold nuclear fragmentation. 

The main common characteristic of cold binary nuclear fragmentation is the 
emission of the final fragments (nuclei) with very low or even zero excitation energy 
and consequently with high total kinetic energy (TKE), which can be provided only 



with compact shapes of fragments at the scission point |T6|, |T^. The cold ternary 
fragmentations should also have compact shapes at scission with deformations close to 
those of their ground states, in order to achieve high TKE values tending to the ternary 
decay energy Qt. It is assumed, that a light nucleus (e.g. an alpha-particle. Be etc.) 
is formed at the surface of the initial nucleus, and that the two heavier fragments are 
formed at the scission configuration in their ground states. 

These cluster like models [|16], [l^ were used successfully to reproduce general 



features of the cold ternary fragmentation. However the scission configuration has 
been built in fact by hands. Therefore it is actual to develop microscopical or semi- 
microscopical approach to this scission-point concept of nuclear fragmentation. 
Recent experimental progress in the production of superheavy elements 



near 



the proposed superheavy "island of stability" |jT9| has stimulated yet the critical 
reexamination of earlier extrapolations pO| , pl| , p2| , p3| of the nuclear shell structure. 
Although this macroscopic-microscopic models quite successfully describe the bulk 
properties of known nuclei, their parameterization needs preconceived knowledge about 
the density distribution and the nuclear potentials which fades away when going to the 
limits of stability. Therefore it is actual to investigate the properties of superheavy 
elements with self-consistent models, like the mean-field models based on the shell 
correction method, self-consistent Skyrme-Hartree-Fock (SHE) WM and relativistic 



mean-field (RME) models (see m, 53, I2§, B^, m and refs. therein) 



There are well developed methods to calculate, in the framework of many-body self- 
consistent approach, static properties of a well isolated nucleus in its ground state. There 
also exists a well developed two-center shell model |^. However, a three-center shell 



model has not been developed yet, except for very early steps [^. Three-center shapes 
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are practically not investigated, in comparison with the two-center ones. There exists 
the generalizations of mean-field models to the case of two-centers [0, but a ternary 
configuration is out of consideration, because of uncertainties to select a peculiar set of 
constraints. 

There also exist a number of calculations for nucleus-nucleus collisions in the frame 
of time-dependent mean-field methods, but an evolution of the cold fragmentation has 
not been investigated yet. 

Therefore, although the principal way to describe nuclear fragmentation in the 
framework of many-body self-consistent approach exists, it is interesting to develop 
other mean-field approaches to analyse these phenomena from different points of view. 

We suggest to use the methods of nonlinear dynamics. These methods gave yet 
the possibility to derive for nuclear physics unexpected collective modes, which can not 
be obtained by traditional methods of perturbation theory near some equilibrium state 



(see [|^, |3^, ^, ^ |33] for the recent refs.). The most importamt reason is that the 
fragmentation and clusterization is a very general phenomenon. There are cluster objects 
in subnuclear and macro physics. Very different theoretical methods were developed in 
these fields. However, there are only few basic physical ideas, and most of the methods 
deal with nonlinear partial differential equations. One of the most important part of 



soliton theory is the inverse scattering method |^, ^, ^ and its applications to the 



integration of nonlinear partial differential equations |^ | . 

In this Letter, we simulate three-center configuration using inverse scattering 
method to solve the single-particle Schroedinger equation, instead of direct solution 
of constrained Hartree-Fock equations. 

2. The framework 

The inverse methods to integrate nonlinear evolution equations are often more effective 
than a direct numerical integration. Let us demonstrate this statement for the following 



simple case. The type of systems under consideration are slabs of nuclear matter ||42 
which are finite in the z coordinate and infinite and homogeneous in two transverse 
directions. The wave function for the slab geomethry is 

1 h'^k'^ 

^kxn(x) = -^=ipn{z) exp(ik_Lr), ek^„ = — -^ + e„, (1) 



where r = (x, y), k_L = {k^-, ky), and Q is the transverse normalization area. 

h' d' 



2 ^ 2^n{z) + U{z)ilJn{z) = e„V^„(z), (2) 

A direct method to solve the single-particle problem (Q) is to assign a functional of 
interaction S (usually an effective density dependent Skyrme force), to derive the ansatz 
for the one-body potential, as the first variation of a functional of interaction in density 
U{z) = U[p{z)] = 6S/6p. Then to solve the Hartree-Fock problem under the set of 
constraints, which define the specifics of the nuclear system. In the simplest case of a 
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ground state, one should conserve the total particle number of nucleons (A), which is 
related to the "thickness" of a slab, via A =^ A = (GApj^/irY^^, which gives the same 
radius for a three-dimensional system and its one-dimensional analogue. As a result, 
one obtains the energies of the single particle states e„, their wave functions ipn{z), the 
density profile p(x) =^ p{z) 

No No 2m 

P(^) = II«nV^n(^)> A=^an, an = —^icp - Cn) , (3) 

n=l n=l ^"' 

and the corresponding single-particle potential. a„ are the occupation numbers, A^o is 
the number of occupied bound orbitals. The Fermy-energy cp controls the conservation 
of the total number of nucleons. The energy (per nucleon) of a system is given by 

_g fl^ ■ ^0 l-OO rinl, ^ No 



it''.iy-t)'-''^^it''ihji:jw)¥^- 



A ' 2m^^^"" ' ' '-' ' "■" ■ " ^""' ■ ' ' -^r.-n— (4) 

Finally, the set of formulas (|l|-§|) completely defines the direct self-consistent problem. 
Following the inverse scattering method, one reduces the main differential Schrodinger 
equation @) to the integral Gel'fand-Levitan-Marchenko equation [^, |39[| 

/•CXD 

K{x,y) + B{x + y)+ B{y + z)K{x,z)dz = 0. (5) 

Jx 

for a function K{x,y). The kernel B is determined by the reflection coefficients 
R{k){ek = h'^k'^/2m), and by the A^ bound state eigenvalues 

N 1 r°° 

B{z) = Y, CliKn) + - / R{k) exp{ikz)dk, Cn = -h'^Kl/2m. 

n=l ^ ■'-°° 

N is the total number of the bound orbitals. The coefficients C„ are uniquely specified 
by the boundary conditions and the symmetry of the problem under consideration. 
The general solution, U{z) = —{h^ /m){dK{z, z)/dz), should naturally contain both, 
contributions due to the continuum of the spectrum and to its discrete part. There 
seems to be no way to obtain the general solution U{z) in a closed form. Eqs. (@),(||) 
have to be solved only numerically. In Ref. [^ we used the approximation of reflectless 
(i?(fc) = 0), symmetrical (f/(— z) = U{z)) potentials. This gave the possibility to derive 
the following set of relations 

U{z) = |^ln(det ||M||) = -— E ^n^li^), 



N 



'4^n{z) = E(M ^)nl\l{z), \n{z) = C„(k„) exp(-K„z), 



n=l 



\niz)\l{z) r^ t ^ /o lA'^^ + '^'lA^/^ 



Mnl[z) = bnl H , Cn[Hn) = (2fi;„|_[J_ M . (6) 

Consequently, in the approximation of reflectless potentials {R{k) = 0), the wave 
functions, the single-particle potential and the density proflles are completely deflned 
by the bound state eigenvalues via formulas (|Tp,(|^). 



Ternary configuration in the framework of inverse mean-field method 
3. Results and Discussion 



In Ref. [Q a series of calculations for the different layers, imitating nuclear systems in 
their ground state was provided. For a direct part of the calculations by the Hartree- 
Fock method, the interaction functional was chosen in the form of effective Skyrme 
forces. The calculated spectrum of bound states was fed into the scheme of the inverse 
scattering method, and the relations were used to recover the wave functions of the 
states, the single-particle potentials, and the corresponding densities. In this note, we 
generalize this method to the case of fragmented configuration, trying to imitate two- 
and three-center nuclear systems. We use here only the inverse mean- field scheme (^. 
The details of the approach and systematic calculations of fragmented nuclear systems 
will be provided in a forthcoming publication. In Fig.l we present the results of the 




Figure 1. The density profiles of tlie light {A k, 20, A ~ 1.0) three-levels 
{N = 3, A^o = 3) model system calculated in the frame of inverse mean field method, 
a) the ground state (solid line); a fragmented two-center configuration (dotted and 
dashed lines); b) a ternary fragmentation of the system into three fragments. 

calculations of three-level (iV = 3, A^o = 3) light {A ?« 20, „4 ~ 1.0) model system 
simulating the ground state (Fig. |l|(a) solid line), and fragmentation of the system into 
two fragments (Fig. |l|(a), the dotted and dashed lines). In the same figure (Fig. |l|(b)) we 
present the fragmentation of the system into three fragments (solid and dotted lines). 
One can see that it is possible to simulate a one-dimensional three-center system via 
inverse scattering method. The following conclusions can be drawn. 

• The density profiles, calculated in the framework of inverse method, are practically 
identical to the results of calculation by SHF method. These results are valid for 
the ground state and for the system in the external potential field. 

• The global properties of single-particle potentials (the depth and an effective radius) 
have been reproduced quite well, but the inverse method yields the quite strongly 
pronounced oscillations of the potential distributions within the inner region, and 
slightly different asymptotic tails of potential. In the framework of inverse scattering 
method, all bound states are taken into account in the calculation of the potential 
(^, but for the density distribution only the occupied states are taken into account 
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(see Eqs. (^). Therefore, the slope of the tails of the potential and of the density 
distributions will we different. 

• We used, the approximation of reflectionless potentials, which gave us the possibility 
to obtain a simple closed set of relations (|), to calculate wave functions, density 
distributions and single particle potentials. The omitted reflection terms {R{k) = 
0)are not important for the evaluation of the density distributions, due to the fact 
that only the deepest occupied states are used to evaluate density distribution (see 
Eq. (^)). The introduction of these reflection terms will lead to a smoothing of the 
oscillations in the inner part of the potential and to a correction of its asymptotic 
behaviour. 

• The presented method gives a tool to simulate the various sets of the static excited 
states of the system. This method could be useful to prepare in a simple way an 
initial state for the dynamical calculations in the frame of mean-field methods. 

4. Conclusions 

Recent experimental progress in the investigation of cold nuclear fragmentation has 
made the development of theoretical many-body methods highly desirable. Modern 
variants of self-consistent Hartree-Fock and relativistic mean-field models give the 
principal way to describe nuclear fragmentation in the framework of many-body self- 
consistent approach. However, the generalization of these approaches to three-center 
case is not provided yet because of existing difficulties to select a suitable set of 
constraints. 

We applied the inverse scattering method to solve the single-particle Schrodinger 
equation, instead of direct solution of constrained self-consistent Hartree-Fock 
equations. It is shown, that it is possible to simulate the one-dimensional three-center 
system via inverse scattering method. 

It is needless to say that the present one- dimensional model is too primitive to 
describe a real three-dimensional fragmentation in nuclear systems. However this model 
may be useful as a guide to understand the general properties of fragmented systems 
and to formulate the suitable set of constraints for the realistic three-dimensional mean 
field calculations of the three-center nuclear system. 



Ternary configuration in the framework of inverse mean-field method 7 

References 

[1] By I.W.Alwarez, as reported by G.Farwell, E. Segre, and C.Wiegand, Phys. Rev. 71 1947 327 

[2] Z. Fraenkel, Phys. Rev 156 (1967) 1283 

[3] Y. Boneh, Z. Fraenkel, and I. Ncbenzahl, Phys. Rev. 156 (1967) 1305. 

[4] A. Sandulescu, A. Florescu, F. Carstoiu, J.H. HamUton, A.V. Ramayya, B.R.S. Babu and W. 

Grciner, Int. J. Mod. Phys. E7 (1998) 625. 
[5] G.M. Ter-Akopian, J.H. Hamilton, Yu.Ts. Oganessian, J. Kormicki, G.S. Popeko, A.V. Daniel, 

A.V. Ramayya, Q. Lu, K. Butler-Moore, W.C. Ma, J.K. Deng, D. Shi, J. Kliman, V. Plhorsky, 

M. Morhac, W. Greiner, A. Sandulescu, J.D. Cole, R. Aryaeinejad, N.R. Johnson, I.Y. Lee, and 

F.K. McGowan, Phys. Rev. Lett. 73 (1994) 1477. 
[6] A.V. Ramayya, J.H. Hamilton, J.K. Hwang, L.K. Peker, J. Kormicki, B.R.S. Babu, T.N. Ginter, 

A. Sandulescu, A. Florescu, F. Carstoiu, W. Greiner, G.M. Ter-Akopian, Yu.Ts. Oganessian, 

A.V. Daniel, W.C. Ma, P.G. Varnette, J.O. Rasmussen, S.J. Asztalos, S.Y. Chu, K.E. Gregorich, 

A.M. MacchiaveUi, R.W. Macleod, J.D. Cole, R. Aryaeinejad, K. Buttler-Moore, M.W. Drigert, 

M.A. Stoyer, L.A. Bernstein, R.W. Lougheed, K.J. Moody, S.G. Prussin, S.J. Zhu, H.C. Griffin, 

and R. Donangclo, Phys. Rev. C57 (1998) 2370 
[7] A.V. Ramayya, J.K. Hwang, J.H. Hamilton, A. Sandulescu, A. Florescu, W. Greiner, G.M. Ter- 
Akopian, Yu.Ts. Oganessian, A.V. Daniel, G.S. Popeko, J.D. Cole, and GANDS95 Collaboration, 

Phys. Rev. Lett. 81 (1998) 947 
[8] E.R. Hulet, J.F. Wild, R.J.Dougan, R.W. Longhead, J.H. Landrum, A.D. Doughan, M. Schadel, 

R.L. Hahn, P.A. Baisden, CM. Henderson, R.J. Dypzyk, K. Summerer, and G.R. Bethunc, 

Phys. Rev. Lett. 56 (1986) 313. 
[9] F.-J. Hambsch, H.-H. Knitter, and C. Budtz-Jorgensen, Nucl. Phys. A554 (1993) 209. 
[10] W. Schwab, H.-G. Clerc, M. Muterer, J.P. Theobald, and H. Faust, Nucl. Phys. A577 (1994) 764. 
[11] Y.X. Dardennc, R. Aryaeiejad, S.J. Asztalos, B.R.S. Babu, K. Butler-Moore, S.Y. Chu, J.D. Cole, 

M.W. Drigert, K.E. Gregorich, J.H. Hamilton, J. Kormicki, I.Y. Lee, R.W. Lougheed, Q.H. 

Lu, W.-C. Ma, M.F. Mohar, K.J. Moody, S.G. Prussin, A.V. Ramayya, J.O. Rasmussen, M.A. 

Stoyer, and J.F. Wild, Phys. Rev. C54 (1996) 258. 
[12] A. Sandulescu, A. Florescu, F. Carstoiu, W. Greiner, J.H. Hamilton, A.V. Ramayya and B.R.S. 

Babu, Phys.Rev. C54 (1996) 258. 
[13] A. Sandulescu, D.N. Poenaru, and W. Greiner, Sov. J. Part. Nucl. 11 (1980) 528. 
[14] P.B. Price, Nucl. Phys. A502 (1989) 41c. 
[15] P. Armbruster, Rep. Prog, Phys. 62 (1999) 465. 

[16] A. Sandulescu, A. Florescu, and W. Greiner, J. Phys. G: Nucl. Part. Phys. 22 (1989) 1815. 
[17] F. Gonnenwcin, and B. Borsig, Nucl.Phys. A 530 (1991) 27. 
[18] In January 1999 element 114 was reported (FLNR, Dubna (the reaction H^Pu + lo^^a). In May 

1999 elements 118 and 116 have been observed in the reaction ^^Kr + ^°*Pb performed at LBNL. 
[19] S. Hofmann, Rep. Prog. Phys. 61 (1998) 639. 
[20] S. G. Nilsson, C. F. Tsang, A. Sobiczewski, Z. Szymanski, S. Wycech, C. Gustafson, I.-L. Lamm, 

P. MoUer, B. Nilsson, Nucl. Phys. A131, 1 (1969). 
[21] U. Mosel, W. Greiner, Z. Phys. 222, 261 (1969). 
[22] Z. Patyk, A. Sobiczewski, Nucl. Phys. A533, 132 (1991). 
[23] P. MoUer, J. R. Nix, J. Phys. G 20, 1681 (1994). 

[24] D. Vautherin, M. Veneroni and D. M. Brink, Phys. Lett. B33, 381 (1970). 
[25] G. A. Lalazissis, M. M. Sharma, P. Ring, and Y. K. Gambir, Nucl. Phys. A608, 202 (1996). 
[26] S. Cwiok, J. Dobaczewski, P.-H. Heenen, P. Magierski, and W. Nazarewicz, Nucl. Phys. A611, 

211 (1996). 
[27] W. Greiner, Int. J. Mod. Phys. E5 (1995) 1. 

[28] M.Bendcr, K.Rutz, P.G.Reinhard, J.A. Maruhn, and W. Greiner, Phys. Rev. C58 (1998) 2126. 
[29] W. Greiner, Invited Talks at Cluster'99, Int. Conf. June 1999, Rab, Croatia, and Int. Conf. 



Ternary configuration in the framework of inverse mean-field method 8 



[30 
[31 

[32; 
[33; 

[34 

[35; 

[36 

[37; 

[38; 

[39 

[40' 
[41 

[42; 
[43 



on Fission and Neutron-rich NucL, June 1999. St. Andrews, Scotland, (to be published in 

Proceedings). 
J. Maruhn, and W. Greiner, Z. Physik 251 (1972) 431. 
J. Hahn, H.- J. Lustig, and W. Greiner, Z. Naturf. 32a (1977) 215. 
J.F. Bergcr,M. Girod, and D. Gogny, Nucl. Phys. A502 (1989) 85c. 
V.G. Kartavenko, Phys. Part. Nucl. 24 (1993) 619. 

V.G. Kartavenko, A. Ludu, A. Sandulescu, and W. Greiner, Int. J. Mod. Phys. E5 (1996) 329. 
V.G. Kartavenko, K.A. Gridnev, J. Maruhn, and W. Greiner, J. Phys. G: Nucl. Part. Phys. 22 

(1996) L19. 
V.G. Kartavenko, K.A. Gridnev and W. Greiner, Int. J. Mod. Phys. E7 (1998) 287. 
V.G. Kartavenko, A. Sandulescu A. and W. Greiner, Int. J. Mod. Phys. E7 (1998) 449. 
I.M. Gel'fand, and B.M. Levitan, Izv. Akad. Nauk SSSR, Ser. Mat. 15 (1951) 309. 
V.A. Marchenko, Dokl. Akad. Nauk SSSR 104 (1955) 695. 
L.D. Faddev, Sov. Phys. Dokl. 3 (1959) 747. 
The Theory of Solitons: The Inverse Scattering Method (cd. by S.P. Novikov, in Russian) (Nauka, 

Moscow) 1970. 
P. Bonche, S. Koonin, and J.W. Negele, Phys. Rev. C13 (1976) 1226. 
V.G. Kartavenko, and P.Madler, Izv. Akad. Nauk SSSR, Ser. Fiz. 51 (1987) 1973. 



